Rydberg molecules in external fields: a semiclassical analysis 
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We undertake a semiclassical analysis of the spectral prop- 
erties (modulations of photoabsorption spectra, energy level 
statistics) of a simple Rydberg molecule in static fields within 
the framework of Closed-Orbit/Periodic-Orbit theories. We 
conclude that in addition to the usual classically allowed or- 
bits one must consider classically forbidden diffractive paths. 
Further, the molecule brings in a new type of 'inelastic' 
diffractive trajectory, different from the usual 'elastic' diffrac- 
tive orbits encountered in previous studies of atomic and anal- 
ogous systems such as billiards with point-scatterers. The rel- 
ative importance of inelastic versus elastic diffraction is quan- 
tified by merging the usual Closed Orbit theory framework 
with molecular quantum defect theory. 
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Gutzwillcr Periodic Orbit theory (POT) represents 
a very successful technique for establishing a correspon- 
dence between quantum properties and classically chaotic 
dynamics. Recently there has been much interest in 
diffractive POT, e.g. j^], for systems with a dynamical 
structure smaller or comparable to one de Broglie wave- 
length. These studies were restricted to a structureless 
point-scatterer, such as a simple s-wave scatterer. It is an 
interesting question, from the point of view of quantum 
chaology, to consider the effect of a diffractive scatterer 
with internal structure, such as for example a two-level 
particle. In fact we show here that this model provides 
us with the key to analyze semiclassically, for the first 
time, the spectra of diatomic Rydberg molecules such 
as H2 in static fields. We find that features such as the 
interference between Rydberg series converging on differ- 
ent energy thresholds depend on a novel type of inelastic 
diffractive trajectory. 

Highly excited hydrogen atoms in static magnetic fields 
represents a well-known paradigm of quantum chaology: 
modulations in the quantum density of states are well 
described by POT However, analysis of the photoab- 
sorption spectra from a low-lying initial state, require 
Closed-orbit theory (COT) jj], a semiclassical descrip- 
tion of a density of states weighted by a dipole matrix 
element. Originally devised for hydrogen, COT was later 
extended to singly excited non-hydrogenic ('Rydberg') 
atoms , accounting for the presence of the ionic core. 
It was found |^ that the core-scattering produces addi- 
tional 'combination' orbits. Non-hydrogenic energy-level 
spectra, even in weak fields where the classical dynamics 
for hydrogen is near-integrable, were found to be signifi- 
cantly more complex, showing short-range level repulsion 
and a proliferation of spectral modulations, a behavior 
termed 'core- induced chaos'. Subsequently, the nearest 
neighbor spacings (NNS) distributions of non-hydrogenic 



atoms in the 'core-induced chaos' regime were shown in 
fact [0 , to correspond to a distribution which was neither 
Wigner (chaotic), nor mixed phase-space, but close to a 
new generic intermediate class known as semi-Poisson 
Q. These findings were investigated experimentally for 
helium atoms in external fields j^]. 

We consider below the form of COT appropriate 
for simple molecules in external fields. The extension 
of COT is obtained by merging the standard closed- 
orbit arguments with molecular quantum defect theory 
(MQDT), the main theoretical approach employed in the 
studies of Rydberg molecules . The success of MQDT 
is based on two ingredients: the separation of the core- 
electron interaction into long-range effects (described by 
Coulomb functions) and short-range effects (described by 
phase-shifts known as quantum defects), and the use of 
different basis sets to account for core-electron couplings. 
Thus, when the electron is near the core, its dynamics 
is coupled to the rotation of the core, the wavefunction 
being described in the coupled frame \a) = \KalaJa)^ 
(A is the projection of the electronic angular momen- 
tum on the molecular axis, I is the orbital momentum of 
the outer electron, quantized along this axis, and J gives 
the total angular momentum). For each a, the quantum 
defect is given by the number and scattering matri- 
ces have a simple representation, e.g. the T matrix el- 
ements are = e^"'^° — 1. When the electron is far 
from the core, the molecular wavefunction is described 
in the uncoupled frame \i) = \Njljmj) {N is the angular 
momentum of the freely rotating core, / is now quantized 
in the laboratory frame along the axis of the magnetic 
field with magnetic quantum number m). These two 
coupling schemes are related by a unitary frame trans- 
formation with elements (j| a). The T matrix elements 
between outer channels j and j'are given in terms of 
by Tjji = {j I a) Ta {a \j') . The frame transformation 
in the field-free case was given by the late U. Fano in a 
classic paper which successfully analyzed experimen- 
tal photoabsorption spectra for H2 from the ground state 
(J = 0, Z = 0) to high Rydberg states {J = 1, I = 1). In 
the presence of a magnetic field of strength 7, the frame 
transformation is modified since J is not conserved, but 
M, its projection on the field axis, is [ |l^ . 

The quantum MQDT treatment, undertaken in an 
inner region where the magnetic field is negligible, is 
matched to the semiclassical one in the following way (see 
Fig. 1). (1) The molecule, initially in its ground state, is 
excited by a laser, giving waves with ^ = 1, J = 1, and 
both A = (phase-shifted by the quantum defect /is) 
and A = 1 (with a phase-shift given by ^u)- (2) As the 
electron leaves the core region, its dynamics gets uncou- 
pled from the core which is left in one of the 4 allowed 
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configurations: A'' = m = and N = 2 m = 0,±1. 
These excited waves are propagated by the full field-free 
quantum molecular Green's function, obtained from the 
Lippmann-Schwinger form of the MQDT wavefunction 

(3) The electron then enters the region in which the 
magnetic field is not 

negligible, where the total wavefunction takes the form 
IV') = En, \NjM - m,) e™^-^Cf ^- Here, V'o":.f ^ is 
similar to the outgoing wave of the electron in COT 0. 
For each core configuration, the electron's wavefunction 
is propagated in the magnetic field region semiclassically 
along classical trajectories k. Two dynamical regimes 
coexist, since the energy partition gives 
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where E^^ is the outer electron's energy and Br the core 
rotational constant (the dependence on m through the 
Zeeman term for the N — 2 levels is very small). (4) 
Eventually some trajectories return to the core. The 
molecular wavefunction is then written as a standard 
MQDT expansion, and the expansion coefficients may 
be obtained by matching the incoming part of the wave- 
function to the semiclassical returning wave. (5) The 
electron is then projected onto the coupled frame, and 
the returning waves interfere with the initially excited 
waves to produce modulations in the absorption spectra. 

These modulations contain contributions from the 
primitive closed orbits k that left and returned to the 
core in one of its 4 possible configurations j and are ob- 
tained by summing terms of the form 
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where Al. and Sj, are the classical amplitude and action 
for trajectory k and energy Ej (Sk includes the associ- 
ated semiclassical phase, containing the Maslov index); 
Ca 1 if A and \/2 if A = 1 is a coefficient stemming 
from the use of the united atom approximation for the 
dipole transitions. (6) The Coulomb field backscatters 
the returning waves which retrace the orbit in reverse, 
giving rise to repetitions without changing the state of 
the core. (7) Core-scattering produces newly outgoing 
waves, which are propagated in the outer region along 
orbits q. These waves now leave the core in state Njirij 
but contain contributions from orbits that had previously 
returned with the core being in state Njimji, with the 
mixing strengths depending both on the classical dynam- 
ics (via the expansion coefficients) and the quantum scat- 
tering process (via the T matrix elements) . These orbits 
return to the core for the second time where they are 
matched once again to a MQDT expansion. The con- 
tribution of these orbits to the oscillatory part of the 
photoabsorption spectrum is of the form 



which clearly singles out the combination of an orbit k at 
an energy Ef with an orbit q having an energy Ej. Fur- 
ther iterations may be performed, but since each process 
of matching the incoming wave to the MQDT expansion 
brings in a factor h}/'^^ fast convergence is expected for 
small h. 

COT for atoms has relied on exact or quasi-exact scal- 
ing properties since the classical dynamics of the excited 
electron depends only on a scaled energy e = E^^'^l'^ . 
The field represents an effective value of Planck's con- 
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Hydrogen in a magnetic field is near- 
integrable for e < —0.5; as the field is increased it makes 
a gradual transition to full chaos at e ~ —0.1. Al- 
though here the dynamics associated with the core in 
states = and 2 cannot be scaled simultaneously, 
scaling techniques may still be applied theoretically to 
unravel the classical dynamics (e.g. [Q). Eq. (0) is 
effectively scaled by giving a 7^/^ dependence to the ro- 
tational constant thereby allowing calculations at fixed 
scaled energies tN=o o,nd tN=2 ■ Full quantal calculations 
of the photoabsorption spectra were obtained by scal- 
ing the method described in [|2 15 1. Note that /-mixing 
brings in (very weakly) the effects of iV = 4,6... states. 
However their contribution is very small and although 
we include these in our quantal calculations, the physics 
remains in effect dominated by the iV = 0, 2 series of 
quantum states. The H2 molecule and its isotopes D2 
and T2 have similar quantum defects. We present the re- 
sults for T2 below since its Br yields the smallest values 
of h in both channels {h - 1/90 and - 1/40 for 
and N — 2 respectively). 

The Fourier transforms of the spectra are shown in Fig. 
2 for M = at a scaled energy ejv=o ~ —0.3, correspond- 
ing to 6^=2 = —2.45, over a range 7^^/^ — 60—120. Over 
this range, Br varies from about twice to half its physical 
value. Each peak is associated with a closed orbit, ei- 
ther a 'geometric' orbit k (i.e. one of the well-known |l^ ] 
closed orbits of the diamagnetic hydrogen atom problem) 
appearing at scaled action Sk or a 'diffractive' combina- 
tion of at least two orbits k and q at actions Sk + S„. 
Whenus = /in, the terms between brackets in Eqs. (g) 
and (H) vanish for Nj = 2 and mj = 0, ±1 (for any j') 
hence only peaks associated with the (eAr=o) closed orbits 
are expected. For zero quantum defects [Fig. 2(a)] there 
is no core scattering (the T matrix vanishes) and only ge- 
ometric orbits are visible. For — nu = 0.5 [Fig. 2(b)], 
elastic core scattering between eN=o orbits produces the 
additional peaks; these combination orbits were seen in 
non-hydrogenic atomic spectra M. (c) shows /is = 0.5, 
/in = —0.5, in which case Eqs. (p|) and ^ predict, as a 
generic feature, the geometric and elastic-scattered peaks 
corresponding to ejv=o to be about 10 times weaker than 
the ones associated with the eN=2 closed orbits. More- 
over, the term given in Eq. (||) vanishes if Nj — and 
Njr = 2,r7ij/ — 0, ±1, totally suppressing inelastic scat- 
tering. Indeed, the regularly spaced peaks in Fig. 2(c), 
are associated with the closed orbits of the near inte- 
grable eN=2 case (the orbits perpendicular and parallel 
to the field for eN=2 — —2.45 have quasi-degenerate ac- 
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tions) and their repetitions. Thus Figs. 2(a)-2(c) display 
simple features characterizing systems in which at most 
only elastic core-scattering is visible in the photoabsorp- 
tion spectrum. 

In contrast, the Fourier spectra shown in Figs. 2(d)- 
(e) contain far more peaks. These arise from inelas- 
tic scattering between the core and the outer electron 
and appear at a scaled action corresponding to the sum 
of an eN=2 orbit (or its repetitions) and an epf^o or- 
bit. The shortest of these combination orbits occurs at 
S = 1.71 and combines the Vi 'balloon orbit' |^ at 
e = —0.3 [peak '1' in Fig. 2(a)] with the primitive or- 
bit 'i' at e = -2.45 [the first peak in Fig. 2(c)]. Note 
that two or more peaks having the same action but aris- 
ing from different combinations will interfere construc- 
tively or destructively depending on the relative phases 
which contain both quantum scattering and classical or- 
bit terms. However generic, peak-independent features 
are governed solely by the quantum defects. For exam- 
ple, when = 0.5, /in — 0, the term between brack- 
ets in Eq. (||) is almost negligible when Nj — Nj' = 0, 
thereby effectively suppressing elastic scattering. This ex- 
plains why in Fig. 2(d) there are no peaks correspond- 
ing to elastic core-scattering of the e7v=o orbits. Fix- 
ing /is = 0.5, the height of these elastic peaks decreases 
monotonically when /in is decreased from 0.5 to 0, as 
expected by our semiclassical formulas. For /is = 0.22, 
/in = —0.06 [Fig. 2(e)], the quantum amplitudes linking 
the different channels have comparable magnitudes, and 
both elastic as well as inelastic combinations are seen. 

We stress that the classical dynamics underlying Figs. 
2(a)-(e) is identical, and the very different plots con- 
firm the importance of the type of core scattering in the 
recurrence spectra. For example the CO molecule has 
equal quantum defects and its recurrence spectra should 
be similar to Fig. 2(b) (elastic scattering only). Fig. 
2(e) on the other hand corresponds to H2 which has the 
well-known quantum defects /is — 0.22, /in — —.06 at 
equilibrium intcrnuclear distance. 

We have also calculated the NNS distributions for dif- 
ferent values of the quantum defects at eat^o = —0.5. 
The results are given on Fig. 3. The distribution for 
A*s = /in = [Fig. 3(a)] appears as a superposition 
of 2 Poisson distributions, the first one concentrated 
at short-spacings (s ~ 0.3), the second one extending 
its tail at longer range. Note that although the Ryd- 
berg series do not interact, the overall spacing distribu- 
tion is not equivalent to a superposition of 2 indepen- 
dent spectra (in which case the resulting distribution is 
well-known |l^). Nevertheless, it was possible to ob- 
tain an acceptable fit to the phcnomenological formula 

P(^) = t(A^)^l(A,-) + f (^)P2(P2,.), With 

X^i Pi — P where p is the total density. Pi gives the 
dominant distribution at short-spacings and vanishes for 
s > 1, whereas P2 gives the distribution which extends 
at larger spacings, with mean spacing p^^ . A is treated 
as a fitting parameter with X > p. Thus for small spac- 
ings, P{s) describes the superposition of Pi and P2 with 



the limiting cases P{s) Pi{p,s) when A — > p, and 
P{s) ~ 6{s — X~^)pi/p + P2{p2, s)p2/ p when A ^ 00. 

Pi{s) was always taken to be a Poisson distribution, 
with A ~ 4. This is the consequence of the weak field 
experienced by the N = 2 levels, which remain clustered 
around their field- free position. However, P2{s) does de- 
pend on the scatterer properties, i.e. on the quantum 
defects. For /is = /in = 0, P2 follows a Poisson statis- 
tics (the tail of the solid line given in Fig. 3(a)). For 
PT, = Pn = 0.5 (b) the tail is in excellent agreement with 
the Wigner surmise, as is also the case for p^ = 0.5, 
pu = —0.5 (c). On the other hand, for /is = 0.5, pn — 
[Fig. 3(d)] and for /is = 0.22, pn = -0.06 (e), P2 is best 
approximated by a semi-Poisson distribution. Globally, 
P2 goes from Poisson to Wigner, through intermediate 
statistics, as the intensity of the scatterer is increased. 
This intensity is given strictly by a multidimensional ma- 
trix obtained from the frame transformation as explained 
above [roughly, the intensity goes from to a maximum 
value as (|/is| , \pn\) goes from (0,0) to (0.5,0.5)]. With 
our simple fit, the specific effects of inelastic scattering 
are not quantified. As known, 2 Rydberg series strongly 
interact when the levels cross each other, and we may 
thus expect the inelastic effects to appear at rather short 
spacings. We do remark however a smaller value of A 
for the cases (d) and (e), where inelastic scattering is 
stronger, which indicates a larger mean level repulsion 
for the distribution modelled by Pi. 

In conclusion, we have shown that the relative contri- 
bution of geometric and diffractive trajectories produced 
by elastic and/or inelastic scattering on the core are es- 
sential to understand the character of recurrence spectra 
of different simple molecules in external fields. The spec- 
tral properties for these molecular systems, containing 
an effective 2-level scatterer, are dominated by 2 spacing 
distributions. 
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FIG. 1. Simplified view of the dynamics of photoabsorption 
of H2 in a magnetic field. The numbered steps are detailed in 
the text, but in brief: (1) A ground-state diatomic molecule 
absorbs a photon. (2) Outgoing waves with 1 = 1 leave the 
core region; the dynamics of the corresponding excited elec- 
tron uncouple from the core leaving it in a superposition of 
rotationally excited {N = 2) or not excited {N = 0) states. 
(3) In the extended region outside the core, the outgoing elec- 
tron dynamics is described using classical trajectories. (4) 
Some trajectories return to the core carrying returning waves. 
(5) The returning waves interfere with the outgoing waves 
producing modulations in the spectra. Returning waves are 
cither (6) back-scattered bj' the Coulomb field ('geometric' 
scattering in the terminology of diffractive POT) producing 
higher repetitions (7), or are core-scattered ('diffractive' scat- 
tering) without exchange of energy with the core (elastically) 
or with exchange (inelastically). The latter produce spectral 
modulations due to combinations of trajectories with different 
scaled energies (e.g. ejv=o and ejv=2) and different classical 
regimes. They are not seen in atomic spectra and are a new 
feature, seen only in the molecular spectra. 

FIG. 2. Fourier transformed quantum spectra for a set of 
diatomic molecules characterized by different quantum de- 
fects ^iB,Mn, showing the relative importance of elastic and 
inelastic diffractive trajectories, (a) fis./in = 0,0: no core 
scattering; peaks correspond to N = series (molecular core 
not rotationally excited / outer electron highly (n ~ 90) ex- 
cited) hence the COs correspond to those of a hydrogen atom 
at e = —0.3. (b) Shows only elastic core-scattering, within the 
N = series, (c) Spectra dominated by COs of N = 2 states 
(molecular core rotationally excited / low excitation {n ~ 40) 
electron) and elastic core scattering combining these trajec- 
tories, (d) Shows trajectories of both N = and N = 2 
states. New core scattered peaks are strongly dominated by 
inelastic diffractive trajectories combining = and N = 2 
trajectories, i.e. the outer electron exchanges energy with the 
core, (e) //E,A'n correspond to H2; shows mixture of elastic 
and inelastic core-scattered orbits. 

FIG. 3. Spacing distributions (NNS) for different values of 
the quantum defects [linear scale for (a), log scale for (b)-(e)]. 
The solid line represents the fit P{s). 
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